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Dark Z/dark photon (Z') is one candidate of dark force carrier, which helps to interpret the 
properties of dark matter (DM). Other than conventional studies of DM including direct detection, 
indirect detection and collider simulation, in this work we take flavor physics as a complementary 
approach to investigate the features of dark matter. We give an exact calculation of the new type of 
penguin diagram induced by Z' which further modifies the well-known X, Y, Z functions in penguin- 
box expansion. The measurement of rare decays B —>■ fi~ and Bs —>■ p~ at LHC, together 

with direct CP violation e'/e in /sT —>■ tttt as well as ^ p'^p~, are used to determine the 

parameter space. The size of coupling constant, however, is found to be 0(1) which is much weaker 
than the known constraints. 
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I. INTRODUCTION 


Dark matter (DM) constitutes about 27% of the energy-matter budget of universe, significantly more than 5% of 
baryonic matter[I|. Nevertheless, the exact nature of DM is kept mysterious so far. Yet it is still unclear whether 
DM can be described by scalar, fermion, vector or even graviton. 

Among various DM candidates, dark photon has been the one of particular interest. The idea was initialed in 1980s 
0, and developed in recent years 0[i|. Suppose there is an extra U{\)£, group besides SM gauge group, under which 
all the SM interactions are invariant. The gauge boson of U{1)d, named dark photon, interacts with SM U{1)y gauge 
boson via a kinetic mixing term. It helps to explain the astroj^ysical observation of positron excesses 0, as well as 
other astrophysical phenomenology such as supernova bounds]^ and Big Bang Nucleosynthesis^. The direct search 
of dark photon, for example XENON 100, has put a very strong constraint 0. For e“''e“ collider, a recent search 
performed at BaBar shows a null result, neither finds nor rules out dark photon 0. A further experiment, the Heavy 
Photon Search (HPS) experiment located at Jefferson Lab [3, is designed to search dark photon in the mass range 
20 MeV to 1 GeV as well as the related coupling. However, it has been realized recently that a simple dark photon is 
not favored by 3.6(T deviation of muon anomalous magnetic moment [HI- Given the fact that dark photon model is 
the extreme case when the parameter to describe Z-Z' mass mixing is closed in a more generic dark Z model, it is 
necessary to extend dark photon to dark Z which is the working frame of this paper. 

Flavor physics is not only taken as a platform for precise test of SM, but also plays an important role in indirect 
search of new physics (NP) beyond Standard Model (SM). Great progresses have already been made since LHG runs. 
For example, it was hoped for decades that NP might exist in zero crossing point of the differential branching ratio 
of H —>■ K*pL'^^~ (see for example [l0) but finally turned out tiny NP effect (l3j| . Taken DM theory as one kind of 
ordinary NP theory, flavor physics then would provide as a complementary way for conventional approaches of DM 
study, including direct detection, indirect detection and collider production. Based on the great achievement made in 
Run I, the LHG has already started its Run H in 2015. Then it would be interesting, and timel y, t o connect together 


li- 


these two different fields, dark matter and flavor physics. A similar effort can also be found in 

The paper is organized as follows. In section 2, we will briefly set up the dark Z model. An exact result of Z' 
penguin and further modifications to X, Y, Z functions are given in Section 3. Some typical processes which are 
affected by the Z' including the recently measured B —^ K*ij,~ and Bg 




as well as e'/e and Kl 


are discussed in section 4, also the relevant formulas are given therein. In section 5 the obtained numerical results are 
shown, based on which we will make a discussion. The details of Z' penguin calculation can be found in appendix. 


II. THE MODEL 


Suppose there exits an extra U{1) group, other than SM t/(l), what it brings in phenomenology is an interesting 
question. It was considered how electromagnetic charge is shifted by this extra C/(l) group in the initial paper 0. 
Until recently it becomes popular to take this U{1) gauge boson as a DM candidate. 

Under the dark group, notated as U(1)d, all the SM interactions are invariant . The connection between dark 
photon with SM particles is from a kinetic mixing term, leading to the effective Lagrangian 0 


where dw is Weinberg angle, Z' and B are dark photon and SM B field with the corresponding field strength 


Buy = duB^ - d^Bu, Z'=duZl-dyZ' 


( 2 ) 


and the mixing of gauge bosons is mimicked by parameter e, which is supposed to be small and need to be determined. 
The convention above, in gauge interaction state, is not diagonalized. By redefining fields as. 



the Lagrangian is then rotated to a diagonal form 

r _ ^ o ^ 


(3) 


(4) 
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Note here the field after rotation with a subscript 0 differs from the one before rotation with a hat. In SM the B field 
can be projected to photon and Z after spontaneous symmetry breaking (SSB). Incorporating Zq, the related neutral 
gauge fields are shifted, 

A = A- eZo 

Zq = Zq + etandv^Zo 

Z' = Z'. (5) 

The rotation does not change the definition of Zq, however, photon and Z field are modified indeed. Due to this 
modification of gauge fields, the interaction between Zq and ordinary matter is induced, which is named as dark 
photon model. 

When Zq-Zq mass mixing is considered, the simple dark photon model is then extended to dark Z model. Generally 
speaking, the mass of Z' could either be added by hand which is called Stiickelberg mechanism[l^ (the origin of 
Stiickelberg photon, for example, is discussed in string theory [T^) or by applying Higgs mechanism, see Q as an 
example. In this paper, we shall adopt the treatment of Z' mass in the latter case, without involving the details of 
the mechanism itself. After the neutral gauge bosons obtain mass after SSB, a further rotation is required after the 
one in eq.([3|) for diagonalising mass matrix 


f Z \ - f cosC -sinC \ / Zq \ 

l^Z'^ - l^sinC cosC JyZ'J 


( 6 ) 


where the rotation angle ( is model dependent and analytically might be complicated, but numerically should be 
small, (for example, see i)- Now combine together the two rotations eq. ([3]) and eq. (|6]) , the modifications to photon 
and Z by dark Z shows 

A^ = A^- eZ; (7a) 

Z^ = cos((Zo_^ — ezZ'^ K, Zo,fj. — £zZ'^. (7b) 

Formally the shift of neutral fields in dark photon model is characterised by two independent parameters e and ez, 
respectively. In fact ez also has a lengthy analytical expression based on detailed model. The equivalent ez defined to 
replace the complicated structure brings the convenience. The dark Z field here and also hereafter is denoted as Z', 
with original SM field denoted with a hat. The two coupling constants together with Z' mass constitute the unique 
3 model parameters of dark Z, which could be measured in experiments. 

Apparently the interactions between Z' and SM particles are simply induced by the shifted neutral gauge field. 
Explicitly, we show how Z' couples to SM fermions. 


^Z'ff 



2 cos 


K 


in which the SM electrical current and weak neutral current are 


( 8 ) 


J^m = Qfirf, 

Jnc = (Tsf - 2Qf sin^ ew)frf - (9) 

where / stands for fermions with corresponding electric charge Qf, isospin T^f = ±i. With both vector coupling and 
axial-vector coupling, Z' behaves as a light version of Z and heavier version of photon. For the coupling of Z' and 
other gauge bosons, it has both “Z component” and “A component” sized by e and ez respectively. 


III. THE Z' EFFECT IN FCNC PROCESSES 

As current energy frontier, the LHC brings plentiful opportunities for flavour physics which dominated by flavour 
changing neutral current (FCNC) processes. In this section, we will investigate these processes of meson physics in 
dark Z model. It is known in SM FCNC processes are induced at loop level. Conventionally the Feynman diagrams 
contributing to FCNC can be classified to three point penguin diagram and four point box diagram. If NP exists, 
the new interaction brought in by NP will modify parts/all of these SM penguin and box contribution. Within the 
Z' model working frame, however, this modification is only applied in photon penguin and Z penguin, keeping box 
diagram contribution unchanged. To make the new effect more distinguishable, we extract these modifications alone 
and name it as Z'-penguin specifically. 
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A. Z' penguin 

We take b —>■ sZ' as an example, noting similar result can be applied to & —^ dZ' and s —>■ dZ' when necessary 
conditions are satisfied. 

In Feynman-t’ Hooft gauge there are totally 10 Feynman diagrams giving contributions to 6 —> sZ'. During our 
realistic calculation, we group two of the external leg corrections and replace them by an effective vertex shown in 
Fig. [TJ with Z' inserted in either two legs. The four self-energy diagrams, effectively two, are then shown as (g) in 
Fig. [21 together with the remaining six ordinary three point diagrams given as (a) to (/). 

G 
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FIG. 1. The effective vertex for & —^ s transition. 
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FIG. 2. The Feynman diagrams contributing to & —^ sZ' in Feynman-t’ Hooft gauge. Fig. (a) to (/) are three point diagrams 
while Fig. (g) is from the correction to external leg where the effective vertex denoted by a cross is explained in Fig[T] 


As mentioned above the full result can be decomposed into “A component” (or e component) and “Z component” 
(or ez component). The e component is same as 6 —>■ sy while the latter one is similar to 6 —>■ sZ. For photon 
penguin, there are two types of effective verteices corresponding to real and virtual photon. In below since we focus 
on semileptonic processes thus only the virtual photon vertex is taken into account. In order to keep the final result 
the same structure as photon penguin, the ez component differs from SM Z penguin by neglecting dipole term 
contribution. An exact calculation, shown in appendix |Al gives the Z' penguin vertex as 

sTf^b\^, = - q^i)(l - -/5)b, ( 10 ) 

2 

in which Ai = VibV*„ (i = u, c, t), Xi = and q is outgoing momentum carried by gauge boson. The vertex function 

Hq^xi) consisting of photon component function Dq^xi) and newly calculated Z component function D^^Xi), are 
characterised by e and ez, giving 


Hq{x) = eDo{x) + ezDo{x) 

4 — I9a;^-b a;^(5a:^ — 2a: — 6) 

D„G) = -- l„i + _ 1)4 


Doix) = -- 


1 


SWCW 


Cw 


34x^ - Ulx^ + U7x - 58 (-3a;^ -b ISa;^ - 27x^ + 19x - 4) In a; 


2I6(-1)3 


36(a: - ly 


-47a;3 -b 237x^ - 312a: -b 104 (3x'‘ - 30a:3 -b 54a:2 - 32x -b 8) Ina: 


108(a:- 1)3 


18(x- 1)4 


(lla) 

(llb) 


(11c) 


with sw = sThOwtCw = cos 9w- Note during evaluating, the light down type quark is supposed to be massless and 
the momentum transfer is small, comparable to light quark mass. We also assume the light Z' mass smaller than the 
mass threshold for muon pair production, which guarantees no Z' resonance is produced when final state of charged 
lepton is muon. Also in this work we will not touch electron and neutrino as the lepton final state for the production 
of Z' when mz' is above electron threshold. Nevertheless, it is safe to neglect mz' in this paper. 
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B. The modification to X,Y,Z function 


In the state-of-the-art effective Hamiltonian approach, physical observables can be factorized into short distance 
(SD) and long distance (LD) contribution. The SD part is treated in perturbative theory while LD hadronic matrix 
resorts to various methods including lattice QCD. The SD contribution gives various combination of penguin diagrams 
and box diagrams and leads to he so-called X, Y, Z functions, which also depends on both the theoretical working 
frame and the calculated physical observable. For example, Y is to characterize the Wilson coefficient in the effective 
Hamiltonian of Bq —>■ and for a more complicated process B —>■ K* pT more functions are involved. In 

the frame of Z' model, we have already discussed a new Z' penguin in previous sector. Before combining detailed 
phenomenology, to include Z' contribution in standard X, Y, Z functions systematically is now our target. 

In SM the X, Y, Z functions are obtained via combing bsj, bsZ, box diagram vertex (see [ 13 ) , and now we need to 
contain bsZ' vertex. The amplitudes of 5 —>• sli (or sb —>■ li) mediated by different gauge bosons are in the form of 


G F Ck — 

iM^ = iX^—= — Do{xi){sb)v-AWv 

V2 27r 


iMz = iXi 


Gf 


iMz' = iXi 


y/2 

Gp a 


2 Co(xi) [vf{sb)v-AiLL)v — af{sb)v-A{LL)A] 


27rs^ 


2Go(x,) 


w 


swHo{xi) , - swHo{xi) , - 


with L = {£, v) and the coupling of Zff and Z' f f 

W = ^3/ - “^QfSw^ «/ = ^3/, 

v’j = e - 2QfSwcw + ^ziTsf — 2QfS^), a'^ = ezT^f- 

For convenience Z and Z' contribution can be put together in a compact form by showing exact final state 

Gp a 




y/2 27rs^ 

Gp a 


•-ZZ' 


’ \f2 27rs^ 

in which we have introduced two parameters 


Go(Xi) [( — 1 -b 4s^)(I + Si)(sb)y-A(^^)v + (1 + S2){sb)v-Ai^^)A'\ 
Go{xi) [(1 -b 52){'sb)v-A{vv)v - (1 + 52){'sb)v-A{vv)A] 


5i = 


awHt){xi) 


4cwGo{xi) [ 1 - 4s 


4civsiv 


+ ez 




S 2 = ez 


swHojxt) 

4cwGo{xi) 


( 12 ) 


(13) 


(14) 


(15) 


Incoperating box diagram contribution, the Wilson coefficients of operators = (s6)y_^ mv, Qfo = {sb)v-AmA 
as well as quark-neutrino operator Q’^_^ = {sb)v-A(vv)v-A can be explicitly extracted as 


GY = A, 


_ -v 

Oio — 


G''’^ = A, 


Gp a 


'J2 27rs^ L 

Gp OL 


s^ • 4 {Zo{xi) + AZo{xi)) - {Yo{xi) + AYvixt)) 


w 


\f2 27rs^ 
Gp OL 

\/2 27rs^ 


Yo{xi) + AYA{xi) 

XQ{xi) + AXo{xi) 


(16) 


w 


in which the X,Y,Z functions are given in the combination of penguin and box diagrams Xq{x) = Gq{x) — 
4Ho(x), lb(a;) = Go{x) — Bq{x),Zq{x) = Go{x) + \Dq{x). The corresponding corrections to X, Y, Z are then given as 


AYv{x) = SiGo{x), 

AYo(a;) = 62 Go(x), 


AZo{x) = SiGo{x), 

AYa{x) = S2Go(x). 


(17a) 

(17b) 


Especially, we find that the modification to Y has two types due to different Dirac structure of lepton pair in the 
operators, while modification to X and Z function are in a fixed way. The modifications can be written in a more 
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explicit form as 


AXo{x) = ^twf-zHo{x) 

XYa{x) = ^twf-zHo{x) 

Hoix) 

Ho{x) (18d) 

In the limit of e —^ 0, (5i = S 2 , the modifications are identical. However, if —>■ 0 (exactly dark photon model case), 

leading to AXq = AYa = 0, then the phenomenology is much more tedious. 


AYy(x) = 
AZo{x) = 




1 - 4s2 


ez 




1 4s^ 


ez- 


tw 

T 

tw 


(18a) 

(18b) 

(18c) 


IV. PHENOMENOLOGY 

The physical observables can be classified into two types in dark Z model. One type relating to box diagram, like 
the mass difference of neutral meson, is not modified by Z'. The other one involving photon and Z penguins, such 
as the direct CP violation in if —>■ tttt, does change. In this section, we will choose several typical processes to see a 
generic effect of Z' on flavour physics. 


A. Bq —^ 


It has been hoped for decades that NP might be unfolded in rare decay Bg —/r+p. . However, no hint of NP 
appeared in Bg —>■ mode from LHC Run I data, given by the full combination results of CMS and LHCb [3 

B{Bg fi+fi-) = (2.8t!(;^) X 10”® (19a) 

B{Bd ^ m+M") = (3.9tl j) X 10-1°. (19b) 

Though Bg —» turns out to be SM-like, there remains a hope for NP in the much rarer mode B 4 —^ (for 

example, see (l^l. 

Due to the precise measurement of decay Bg —>■ is now realistic, one should consider the effect of sizeable 

width difference AP, in Bg-Bg oscillation. The theoretical formula has to be corrected to compare with measured 
branching ratio [2^ which is denoted with a bar, 


B{B° £+£-) 


l-?/s 


l + -4^r y^- 


B{B° £+£-) 


( 20 ) 


where yg = 


AFs 

2rs 


_pC^s) 


2 Ts 


Ae+e- 

^Ar 


ig known = 1 in SM [H, thus 


R' 




+R 




^ ^ ^ (21) 

the latest estimation of parameter yg is yg = 0.069±0.006 given in [2^. Note in the dark photon model, the relation of 
eq. dm does not change. While for the rarer Bd ^ £^£ decay, the effect from oscillation in B^-K can be neglected 
thus we do not take this correction. 

The (uncorrected) SM branching ratio of Bq —7> £^£~ is induced by Z penguin and hence depends on Qio, (see ref. 
[H). Now incorporating Z'-penguin, which gives a similar component as Z, leads to 


B{Bq 


£+£-) = T{Bq)^ 
TT 


a 


\47rs^ 


flmjniBA A - 


Am 


I 


»?eff 


Xt(Yo{xt) + AYA{xt)^ 


( 22 ) 


with ?7efi = 0.9882 ± 0.0024 which takes into account NNLO QCD correction and NLO electroweak correction [^ . 
Apparently the reason why only AYa contributes, is exact with the same reason why photon penguin contribution 
vanishes. 
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B. B ^ fi+fj,- 

The quest for NP in _B —>■ H has been performed for a long time. In the beginning the zero crossing-point 

Qq is of the first priority, however, turns out to be compatible with SM prediction finally from the released LHC 
data. The remaining possibility for NP in this mode, the problem, requires more data to confirm. Meanwhile for 
the B —>■ channel, there is a NP hint, so called Rk problem, which violates lepton universality. 

The theoretical study for this channel has been developed for around 30 years, including multi-loop calculation of 
Wilson coefficients at high energy. The most matured theoretical treatment in low energy to the semileptonic decays 
are based on QCDF. It is not necessary to repeat the whole long story in this work. Instead, we would like to simply 
focus on Cg and Cio to see how the data constrain the NP parameter space. In dark Z model, the modification of Cg 
and Cio due to the dark Z contribution is 

ACg{x,) = 4AZo{x,) - (23a) 

Sw 

AC'io(cr,) = (23b) 

Sw 

Driven by new data many efforts have been made to extract the information hidden inside the two coefficients. In 
a model-independent fit was taken based on B ^ Xs£'^i~, B —?> Xg^, B —>■ K*"f and B —?> we will adopt 

their constraints on ACg and ACio at 2 a Q 

- 1.5 < Re (ACg) < 1.2, -2.8 < Im (ACg) < 2.8 (24a) 

- 1 < Re (ACio) < 1-5, -3 < Im (ACig) < 3. (24b) 

In our scenario, the NP only exists in the change of real part of Cg/iQ. 


C. jKl —^ 


The branching ratio of t /r/I contains LD and SD contribution. The calculation of LD contribution remains a 
challenge in theory. Combining a latest theoretical LD estimation and experimental bound [1^ 0, the constraint to 
SD is 


B{Kl ^ < 2.5 X 10 


-9 


(25) 


The branching ratio from SD (see ref. IH) 


B{Kl )sD = 


Re(Ae) 

IK. I 


Pc{Yk) + 


Re(At) 


1 2 


wYoixt) 


with Xg = V*gVgd{q = C,t) and = (2.009 ± 0.017) x 10“® Pc(Kf) = (0.115 ± 0.018) [13, QCD 

correction factor t/y = 1.012 [^, now is modified as 


B{Kl —t )sD = 


by including the dark Z contribution. 


PcIYk) + ^y^l^ wiYoixt) + AYAixt)j 


IP 


IK 


(26) 


D. 


e' 

e 


Historically two approaches, operator production expansion (OPE) method and penguin-box expansion (PBE) 
method, are adopted for the study of direct CP violation in RT —>■ tttt, which involves all the QCD penguin and 


^ Later in another independent analysis [^ . the global fit combining B —>■ , B —> K*pApi~, B —t KpX'p,~ and Bs —t pYpL~ 

together obtained a similar constraint on Ciq. 

^ The SM prediction is then B{Kl — t P'^P~)so = (0.79 ± 0.12) X 10“®, the experimental value in PDG is 
B{Kl At+M“)exp. = (6.84 ±0.11) X 10-9. 
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FIG. 3. The allowed e — ez parameter space by various experiments if ruz' is less than 2m^ and hence ignored. The meaning 
of colours are given as: grey stands for Kl —>■ orange represents Bs —>■ green denotes ACio, blue stands for ACg 

while the yellow is one example of e'/e with {Re, Rs) ~ (0.7, 1.1) and the black dot remarkes SM case. 


electroweak penguin in SM. For the phenomenology study here, we make use of the simple analytical formula based 
on the PBE methodjs^. By modifying relevant parts due to the dark Z effect, an updated formulcd to depicted e'/e 
is 

6 ^ 

Re— = alm(At) • F{xt) (27) 

e 

where F{x) is given by 


F{x) =Po + Px [Xo{x) + AX (a:)] + Py [Ib(^) + XY{x)] + Pz [Zo{x) + AZ] + PeEo(x) (28) 

and the factor a = 0.92±0.03 [s^, which takes into account the correction due to AI = transitions [sj. Note the 
dark Z modifies most parts of SM F function but keep the gluon penguin vertex F{x) unchanged. The coefficients 
Pi {i = 0,X,Y, Z,E) are given in terms of Rq, Rs 


-JO) 


Pi = r, 




(29) 


We adopt their numerical values for as{Mz) = 0.1185 given in Table 1 of Ref. 33|. For the nonperturbative 
parameters,we adopt the value Rs = 0.6, Rq = 1.1. The former one is obtained from lattice [s^, with the translation 
by Ref. [s^. But a reliable lattice result for i?g is still lack, here we choose 10% deviation from large N result. The 
experimental value with Icr error for s'/e is 


taken form PDG (35j| . 


j ~ Re = (1.66 ± 0.23) x 10■^ 


(30) 


V. RESULTS AND DISCUSSION 

There are limited three free parameters in dark Z model, e and ez- In our working scenario, the light Z' mass 
is ignored thus it might be promising to determine the remaining two by above observables. We take the global fit 


^ We should keep in mind that the original formula was obtained by comparing with two methods in SM Is^I . A more serious formula in 
dark Z model should be given by repeating this work similarly due to different types of AY, which is beyond the scope of current work. 
In the numerically study below, we will take AIV an example. However, we will understand the exact form of AY should keep the 
paper’s conclusion. 
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of Wolfenstein parametrization of CKM matrix as input since Z' does not change SM flavour structure. The other 
related input parameters have been given herebefore. Combining Bg —>■ —>■ and ACg, ACio from a 

global fit of i? —>■ and so on, we plot allowed parameter space in (e, e^) plane shown in FiglSl 

The ranges of e and ez are both shown in the section (—10,10). Due to its large uncertainty SD of —>■ 

as presented in grey region, gives a pretty wide band with | slope. The green band with the same slope and less 
width of grey one, together with the gap out of this band surrounded by a set of hyperbolic curve, is resulted from 
ACiq. Another parameter obtained from B —>■ samely, ACg, constrains parameters in blue colour, which 

is restricted in two sets of hyperbolic curves and extended to slope | direction as well as near vertical direction. 
Apparently the latter part is excluded by ^ (or Cig) while the former one is embedding in green area. 

Remarking in orange, the important Bs —>■ also embeds in ACio but somehow has an overlap part with ACg. 

In principle, B^ —^ can also be included. Considering its uncertainty is larger than Bg —>■ the allowed 

region is then also wider than current orange area thus we do not show it. Here we also add in a typical constraint 
from e'/e by fixing its non-perturbative parameter Rq and i?8- Showing as two narrow hyperbolic curves in yellow, 
it looks e'/e could give a very strict constraint. However, due to the large uncertainty of the two non-perturbative 
parameters, especially Rg; one cannot take the direct CP violation of kaon too seriously. Nevertheless, once Rq is 
fixed within certain precisement in the future, e' je can be taken as an important discrimination to further constrain 
the parameter. 

Aiming at the determination of parameter space of dark Z model, we reach the allowed (e, ez) in the narrow linear 
region 


3 




(31) 


and especially —2 < e < 2, —1 < < 1 are favored. In other words, the order of mixing parameters is constrained 

to be 0(1). Generally speaking, more FCNC processes can be considered to determine the bound, but we can believe 
0(1) should be the typical value from flavour physics. 

Many other works have already put constraints on dark photon model. Though it is not exact same model as our 
working scenario, the obtained mixing angle in dark photon still enlightens parameters in dark Z model. For example, 
based on supernova 1987A the limit on mixing angle could be 0(10“^^ ®) for mz' > 2me [1]. A recent work from 
direct detection experiments such as XENONIO and XENONIOO, for the absence of an ionization signal, puts a even 
more stringent limit on e down to 0(10“^®). One may expect if ez is added, the bound of e will also be changed from 
astrophysical observables and direct detection experiements. Meanwhile, the working scenario in flavour physics can 
also be modified to allow Z' resonance production Q, leading to a possible different constraint on mixing parameters. 
In any case, given the bound from flavour physics in current scenario 0{1), one may not expect a dramatical change 
(a more than 10 orders of magnitude) happens unnaturally. 

In the existence of dark Z, we have investigated its effect in flavour physics especially by connecting the effect to 
the newly measured processes at LHC. However, the obtained bound, 0(1), may not compete with the corresponding 
one in traditional dark matter study. 
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Appendix A: The calculation of Z' penguin 

The effective bsZ' vertex is in the form of 


r^(fci, fcg) = eAi ——- g^7^)PL -b Fg , (Al) 

z(47rj rriyy 

which satisfies Ward Identity approximately in the limit of mz' —>■ 0. In the phenomenology study of this work, we 
only focus on Fi contribution, which has the relation with iLg in section IHI Al 


Ho(x) = -2Fi{x). 


(A2) 
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For convenience, in following calculation the effective vertex (EU can be rewritten as 


r'^(A:i,fc 2 ) =eA. 


92 


2(47r)2m^ 


[clT^ + + CsT^] 


(A3) 


by introducing Cf, = Ai + F 2 ,Cs = F 2 — Fi,cl = —F 2 and Particularly we 

have 


Ai = i(cb-Cs). (A4) 

In the calculation, the Fi function may be projected to 

Fi = eFil^ + ezFx\z (A5) 

Following we will evaluate Ff component by component. 


1. The e component 


The contribution of e component is exactly same as & —> 57 case. Here we recalculate the photon penguin contribu¬ 
tion. In Feynman-t’ Hooft gauge, with the approximation iris = 0,q^ = 0 and making use of quark on-shell condition 
as well as unitary CKM relation, we have the total contribution as 






UK 


(A 6 ) 


with the individual amplitude from each diagram 


Ma\=-, 


3 L 


(—2Coo + rri^Co — TOj(Co -I- 2Ci -I- C 2 -f Cn -I- C 12 )) 7l 


+ 2(Cii + Ci)Tjj + 2(Co + Cl -f C 2 + Ci2)T's 


Mb =- 


Me\=- 




3 'm-w 


—2Coo 2 C'o ~ Wfc (Cii + Cl + C 12 -l- Co) ) Tl 


Md =- 


+ 2C11T6 + 2 (Ci 2 -C 2 )r« 

(12C00 + m^( 3 Ci -|- 2C2 -l- 2 Cii -|- 2C12)) II 

+ (4Cii + 2 Ci)T6 + ( 4 C 21 - 2Ci - 2C2)Ts 

2 CooTl + ( 2 Cii -|- 3 Ci+Co)Tb + (2C21 -|- C1-I-2C2 -I- Cq)Ts 


TO? ■■ 


Me\ = 


ruw 

m^tCon 


Mf\= (TO?Co+TO?Ci)Ti+ 2 C 2 r, 


Mo = 


(Ho + Hi) 2 + 


mr 




mr 






Tl 


(A7) 


in which Bi,Cj (* = OCc = 0,1,2,11,12,21,22) are Pasarrino-Veltman integralsjs^, and the positions for vari¬ 
ables are assigned as B{ml,m‘^,mf), Ca,f,(TO?, 0 , 0 , to?j^,to?, to?) for a,b case, and Cc,d,e,f{'ml,0,0,mf,m'^,m‘^) for 
c, d, e, / case. Note in the calculation, the light down type quark contribution in Goldstone-quark-quark vertex cannot 
be neglected. 

The PV functions can be reduced into basic scalar function Ho, Co, and we further perform Taylar expansion up to 
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mf. Sum up all the contribution together, we have 


92 


'7 *2(47r)2m^ [V 3 


1 


Tl 


~77 ~ 2Ag ) — y + Ci,T]^ + CfeTfe + CgTs 






X( (2 — 3xt) Inxt 22X( — 153X( + 159xt — 46 
2(xt - 1)4 36(xt - 1)3 

(3xt — 3xt + 36X( — 32xt + 8 ) Inx* 19xj — 222x2 + 165xt — 34 


18(xt - ly 


108(xt - 1)3 


- 


(—3X( + 57x3 _ Y2xf + 32xt — 8) Inxt 113x3 — 696x3 _|_ _ 242 


18(xt - 1)4 


108(xt - 1)3 


The first term in eq. (IA 8 al) gives zero contribution after applying the unitary triangle relation ^ 
Also we have 

I I (3x4 — 30x3 + 54x3 — 32xt + 8 ) Inxt —47x3 + 237x3 — 312xt + 104 
- 9(xt -1)4 54(xt - 1)3 

= “2 ^Zlo(2;t) + —^ 


Fi|^ = -Do{xi) 

in which DQ{xi) is vertex function of virtual photon in 6 —>■ 57 . Note in my calculation, the Dq^Dq 
one in [13, up to a minus sign, for the different convention in QED vertex. 


2. The ez component 


The total contribution for b —>■ sZ is 

*^1 z = * E 7^ U 

2(47r)2 cw 3^ 

3 

Compared with b ^ case, we only need to recalculate Fig.(a) and Fig.(b), 

2 


^a\z ~ ^6(^11 + C 12 + 2Ci + C 2 + Cq)]T^ 

+ 2(Cii + Ci)Ty + 2 (Ci2 + C 2 + Cl + Co)ri^} 

- i • 2{[-2Coo - mliCn + C 12 + 2Ci + C 2 + Co)]7 

+ 2(Cii + Ci)r,'^ + 2 (Ci2 + C 2 + Cl + Co)ri^} 

{[- 2 C 00 + i + m^Co - m3(Cii + C 12 + Cl + Co)]7 
+ 2CiiT,^ + 2 (Ci2 - C 2 )Ty} 

{[m?Co - m3(Co + Ci)]7 - 2C2ri^} 


I _2 2 m3 


1 . 

2 m^ 


while the other contributions are obtained by the following replacement 

Fig.(c) : e g 2 Cw 


Fig.(d) : e^g2 


1 - 2s3 


w 


2civ 


Fig-(e),(/): e^- 


92Sw 
C w 
1 1 


’IT 


(A 8 a) 

(A 8 b) 

(A 8 c) 

(A 8 d) 
,i • constant = 0 . 

(A9) 

(AlO) 
differs from the 

(All) 


(Al2a) 

(A12b) 


Fig.(5) : -Qb 


(A13) 
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with the same PV function convention as in 6 57 case. 

In this work, since we only consider the (< 7 ^^ — 9 ^ 7 /^) term, thus only the coefficients of and are of interest. 
Especially they are listed individually, 


‘■w 


pM _ 

a z 


—5x^ + 22xt — 5 (1 — 3xt) In a;* 


18(xt-l)3 ^ 3(a;*-l)4 ^ 

20a;( — 7xt — 7 (fix^ — 6xt + 1) Inxt 


2(5x| — 22xt + 5) 4(3a;t — 1) In xt \ 


27(xt-l)3 + 9(a:i-l)4 


18(xt - 1)3 3(a:t - 1)4 

2 ^ 2(20x4 ~ ~ 4(6x4 ~ 6x4 + l)lnx4^ 


27(xt - 1)3 ' 9 (x4 - 1)4 J 

2 T.MI _2 f 2(2x^ - 7X4 + 11) , 4X4 lnx4 ^ 

^wh\z 27(1-X4)3 ^9(x4-l)4j ^ 

Xtixt-5) Xtinxt 


rpfl 


+ 


4(xt - 1)^ 2(xt - 1)3 

2 j X4(11 x 4 — 43x4 + 38) 2(3x4 — 4) lnx 4 ^ 

+Sr 


mw^c\z = 


_ 

27(x4 - 1)3 9(x4 - 1)4 )' “ 

— 17 x? — 8x4 + 1 X4(x4 + 3)lnx4 2 / 17X4 + 8x4 — 1 X4 (x4 + 3 ) lnx4\ 


3(x4 - 1)4 + 


18(x 4 - 1)3 
—38x4 + 43x4 — 11 2 x|(4x 4 — 3) lnx 4 


9(x4 - 1)3 


3(x4 - 1)4 


2 I 2 x^(4x4 -3)lnxt , 38x^ - 43 x 4 + 11 ^] ^ 4 . 
3 (x 4 - 1)4 + 9 (x 4 - 1)3 


Tt 


rpfl 

j\ ^ 


''Xl-W^dlz ~ 


xt{x‘l — 8 x 4 — 17) x^{xl — 3 x 4 + 6 ) lnx 4 


72(x4 - 1)3 


2 {xtixl - Sxt - 17) 0,4 

-\-Sw 


12(x4 - 1)4 

xi (x| — 3 x 4 + 6 ) In X 4 \ 


36(1 -X4)3 6(x4 - 1)4 

X4(23x 4 — 22x4 — 13) (X 4 + 3x4 — 6 ) lnx 4 


72(x4 - 1)3 


12(x4 - 1)4 


) I n‘ 


2 /x 4 (— 23 x 4 + 22 x 4 + 13) xf (xj + 3 x 4 — 6 ) lnx 4 \ 

' 36(x4 - 1)3 6 (x4 - 1)4 ) 


1 - 3X4 


'^W^e\z 

™'^^/I^="'^(-2(x4-1)2 (x4-1)3 

'^W^g\z 


a;4lnx4 1 ^ 4 , 


rpg, 


Sum up together, we have 


Cb\z 


—35x3 _|_ x 563;2 _ 213x4 + 20 (—3x4 _|_ g^s _ -|- 28 x 4 — 4) lnx4 

216(x4 - 1)3 36(x4 - 1)4 

2 f 19 x 4 — 222 x 4 + 165x4 — 34 ( 3 X 4 — 3 X 4 + 36 x 4 — 32x4 + 8 ) lnx 4 

V 108(x4 - 1)3 ^ 18(x4 - 1)4 


(A16) 


Cs 


z 


-103x3 + 438x2- 507x 4 + 136 , (3xf - 27x? + 36x? - 10x4 + 4) lnx 4 
216(x 4 - 1)3 36(x4 - 1)4 

2 /113xf - 696x? + 789x4 - 242 (-3x^ + 57x? - 72x| + 32x4 - 8) lnx 4 

V 108(xi -1)3 + 18 (x4 - 1)4 


(A17) 
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Now we obtain the Dq function as 




34X( — 141a:( + 147xt ~ 58 ^ ~ + 19xt — 4) Inxt 




—ilxf + 237Xf — 312a::t + 104 ^ (3xj — 30Xf + 54xf — 32xt + 8) lna;t 




(A18) 


Take into account Dq and Dq together, the vertex function for Z' penguin Hq is produced. 
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